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Introduction 


This report covers work done at LARS/Purdue under the Landsat-D(*0 
Image Data Quality Analysis Contract for the period August 9 through 
November 9, 1982. In this period, four-band Thematic Mapper data vere 
received on September 9 over the Detroit area and seven-band TM data 
from northeastern Arkansas were received on September 29. Reformatting 
software work Md begun early so that the data could be viewed and pro- 
cessed at LARS/Purdue shortly after receipt. Tn general, the TM data 
appeared to be of very high quality with the expected high resolution 
and very low geometric and radiometric error content. No Laud sat-4 MSS 
data were obtained in the period. Numerous readily applied analysis 
procedures were carried out to provide a "quick look" analysis to iden- 
tify obvious problems prior to detailed parameter evaluation. These 
results are presented here. 


Problems 


The most significant problem is the lack of an acquisition of a 
TM/MSS frame over the prime test site, Chicago, to date. An alternate 
prime site is being considered for initial Investigations. This would 
be our Priority 4 site including Jacksonville, FL, the Atlantic Coast, 
and numerous interstates. An acquisition is available from October. 


Significant Results 

This section contains the "quick look" analysis for the Detroit 
and Arkansas TM data. The work is divided into geometric and radiome- 
tric considerations. In general, the TM data appeared visually to be 
of very high quality; however, several errors were noted in some bands 
and these are briefly identified. 


Geometric Evaluation 


In the Detroit scene, several geometric problems were noted in 
Band 2. Significant scan line misregistration was observed in this 
band on a nominal cycle of 16 lines. A line correlation algorithm was 
developed by LARS for use in this study and was applied to the data to 
analyze channel-to -channel and band-to-band registration. This algor- 


-2 


ithm was appllad to a blook of data near downtown Detroit. Tha 
channel -to-channel and band -to -band results are all very good, averag- 
ing .02 pixel except for one deteotor of Band Two.* This has an indi- 
cated Misregistration of 1.1 pixels to the east. In addition, whan 
Band 2 was blink coopered to Bands 1 and 3, blocks of imagery in Band 2 
of approximately 16 lines by 128 columns appeared to be vertically mis- 
registered by one to two lines. Other misregistration shifts of Band 2 
relative to others were visually observed but not analyzed due to 
reoelpt of the seven-band Arkansas test scene. 

The problems noted in Band 2 of the Detroit scene were not 
bserved in the Arkansas data and it was assumed that these were obvi- 
ous preprocessing problems which had been corrected and were not faults 
in the sensor Itself. A block of 512 by 512 points in the Arkansas 
data was defined as a test site. This block surrounds NASA AgRISTARS 
test site Segment 306 which we have seme familiarity with. Detailed 
line correlations were run on the bands and channels of this Arkansas 
data subscene and the results are in Table 1 . The rows in the table 
are correlation results for each detector for the reflective bands for 
correlation over 100 pixels across a scan line and the result averaged 
down 32 lines to encompass 512 lines. The columns of the table are 
correlations between bands as noted at the top of the column. The bot- 
tom two rows are averages and standard deviations over all deteotors 
for uach band pair. Each detector now consists of a row of mean values 
and a row of standard deviations. 

The results for the within spectral region bands (i.e. , 1,2,3) or 
4,5,6 are extremely good. Averages as low as .01 pixel are seen. For 
correlations between band groups, the averages Jump to the .5 range. 
All the bands indicate a significant misregistration of the thermal 
band with averages between 1.52 and 2.17 being observed. Individual 
detector shifts of over three pixels were noted. Two questions arise 
when carrying out these correlations : One is that as correlation bet- 

ween two images decreases, the variance or the misregistration estimate 
increases; so some of the Indicated misregistration is due do this 
effect. This problem will be better quantified in future reports. The 
second • uestlon is the effect of the cubic convolution interpolation 
resam; x+ng on registration estimates. The geometric correction process 
uses this resampling method and it may blur registration differences 
between detectors. The "A" tape data will be free from this effect. 
Block correlations for the entire 512 by 512 area (Table 2) were com- 
puted to aid in determining the likely effect of misregistration esti- 
mation variance on the values in Table 1 . It is interesting to note 
that the visible bands are rather highly correlated with the third 
reflective IR band and to some extent the second. Correlation is low 


• We (LARS/Purdue) have interchanged 
tapes so that the thermal is Band 7 to 
cally increasing. This should be noted 


the last two TM bands on our 
make the wavelengths monotonl- 
ln all future disousslons. 


ORIGINAL • 

Of POO* QUALITY Cabh L 


Means and Standard Deviations for Arkansas TM Data. 
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Table 2. Block Correlation Tor 51 2 x 512 Points In TM Data 

for Arkansas Sits. 


Band 

1 

2 

3 

4 

5 

6 7 

i 

1.0 






2 

.92 

1.0 





3 

.95 

.94 

1.0 




4 

-.50 

-.46 

-.58 

1.0 



5 

.53 

.54 

.53 

.14 

1.0 


6 

• 

CD 

O 

.79 

.84 

-.37 

.81 

1.0 

7 

OO 

on 

• 

.36 

.41 

-.52 

.19 

.43 1.0 


with the thermal IR. The low correlation could be due to the misregis- 
tration of the thermal band, although the effeot may be small because 
of the 120 meter resolution. 

Band pairs having the approximately .5 misregistration are (1.5), 
(1,6), (2,5), (2,6), (3 >5), (.>,6), (4,6). This pattern definitely sug- 
gests that Bands 5 and 6 (1.55 - 1.75 and 2.08 - 2.35 ym) are unregis- 
tered with respect to 1,2, 3, 4 by about half a pixel. The correlation 
coefficient is quite high between the visible bands and the 2.08 - 2.35 
m band, so the misregistration estimate should be reliable. LARS 
Bands 5 and 6 themselves appeared to be well registered. 

Visual blink comparisons were run on the thermal band using the 
LARS Comtal digital display and it was concluded that the misregistra- 
tion was four lines and columns (right and down). The line misregis- 
trations could account for the across-track numerical correlation not 
producing a value of four columns of misregistration. We feel that 
this is a software error related to the geometric correction and cubic 
convolution and should pose no big problem to oorrect. LARS is regis- 
tering the thermal band for the test site based on these results. 

Considerations were made of the bow-tie and scan-line nonlinearity 
effects but no results were obtained on these questions in the quarter. 


Radiometric Evaluations 


Visual and numerical analysis was carried out on the Detroit TM 
data. Band 2 demonstrated radiometric errors In the form of dark mul- 
tiplxel spots randomly distributed over the frame. The dots have a 
site of up to 4 by 4 pixels and values of 10 to 40. 

Receipt of the seven -band Arkansas frame prompted a more detailed 
analysis of the radiometric characteristics. Visual Inspection was 
first carried out on all the bands and It was noted that the bad data 
spots were not present In Band 2 and all the bands appeared to be of 
very high quality, except for the thermal which had dot problems and 
rather low contrast. 

The first numerical investigation was computation of means, vari- 
ances, and histograms for a block of data surrounding Segment 306, as 
discussed above. The intention was to evaluate the channel-to-channel 
calibration of the 16 detectors for the first 6 bands and the 4 detec- 
tors for the thermal band. He note again that the cubic convolution 
may tend to alter results as compared to what would be seen in the data 
from "A" tapes. Histogram plots for one detector are presented In Fig- 
ure 1. The data ranges for the visible bands are small. The first IR 
band is the only one which uses most of the 256 bln range. The thermal 
band has a limited range with a standard deviation of only 5.6. The 
data distributions appear to be ske ed to the higher values for all 
bands. This was noted in most aircraft and Landsat 1,2,3 MSS data In 
past investigations. Explanation for this will not be attempted here. 

The means and standard deviations for lines taken at 16-line 
intervals were computed in an attempt to estimate channel-to-channel 
calibration. These results are presented in Table 3, a through g. The 
means were all within .2 count of the grand mean for all detectors of 
all bands, except for four detectors of Band 4 which deviated from the 
grand mean by up to .49 count. The specification on channel-to-channel 
calibration Is that deviation be less than the RMS noise level divided 
by 4. The radiometric accuracy specifications for each band are given 
in Table 4. The deviations in the means in Table 3 are due to sensor 
and calibration variations plus variations due to sampling of different 
areas of the scene. The 16-llne interval Is repeated 32 times in the 
512-llne block. Thus there are 32 x 512 s 16,384 samples In each mean 
calculation. In Table 4 we have computed a working estimate of chan- 
nel -to -channel mean deviation tolerance for this quick-look analysis. 
Based on these rather gross assumptions, we conclude that all channels 
of all bands are generally within tolerances, except for about half of 
the deteotors In Band 4. Some, such as detectors 1,2,3 and 11, are 
three times the tolerance. He cannot say for certain whether there Is 
a problem here at this time In bias settings, but there is some Indica- 
tion from these results that there may be. No evaluation was made of 
the variance differences. 
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Figure 1. Histograms for TM bands for one detector for 512 x 512 point area 

from Arkansas test frame. 
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Figure 1 


(Continued) 
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Other rad lone trie observations were made in the therm 1 data. 
Thera are vertical rows of bright dots in this band. The dots have a 
vertical Interval of 17 to 18 lines and are spaced 128 pixels apart 
horizontally. The value of the dots is approx lea tely 157. The size of 
the dots is one pixel wide and up to 3 pixels vertically. 

This radiometric and geometeric quick-look analysis hopefully 
identifies obvious problems and to some extent evaluates the key per- 
formance parameters of registration and channel-to-channel calibration. 


Publications 


No publications were produced in the quarter. 


Recommendations 


Two technical recommendations arise out of this early analysis. 
One is to further check and adjust the registration of Band'i 5, 6, and 
7 if they prove to be out of registration as is suggested by our 
results. The second is to examine the thermal IR processing stream to 
find and correct the cause of the bright pixels in the data. A third 
recommendation is to Interchange TM Bands 6 and 7 so that wavelength 
Increases monotonlcally with band number. 


Funds Expended 

The funds expended on the project are reported perloolcally by the 
Purdue Office of Contract and Grant Business Affairs to the sponsor on 
NASA Form 533M. Tht-se are Issued monthly. Specific disclosure of 
funds expended in this form of technical report is not permissible. If 
a quarterly summary is required, it can be prepared as a separate docu- 
ment by Purdue OCGBA. 


Data Utility 

Data utility conclusions as a result of these initial investiga- 
tions are, first, that Band 2 of the Detroit data has some problems and 
should be used with caution. Similarly, the thermal IR band of the 
Arkansas frame appears to be significantly out of registration and 
should be used with caution in conjunction with the other bands. 

Finally, it was of Interest to examine the dimensionality of the 
new seven -band TM data relative to what had been available from the 


-9- 


Landsat MSS for the past ten years. Principal components analysis was 
conducted of the 512 by 512 block in the Arkansas data using three sets 
of bandst The first four, the first six, and then all seven. We felt 
the first four would reasonably represent the MSS bands as we did not 
have MSS data for this scene. Again we emphasize that our Band 7 is 
the thermal IR. The eigenvectors and eigenvalues for the three cases 
are listed in Table A. The results for the four-band ca ' b are similar 
to those for previous MSS bands even though there Is only one IR bai.d 
and a blue band is present here. There typically are two components 
for vegetated scenes in this wavelength range. Adding the two middle 
IR bands Increases the dimensionality to three, indicating significant 
correlation of these two bands } but as a pair, they are relatively 
uncorrelated with the visible and n ear IR regions. Adding the thermal 
Increases the dimensionality to about four. These are expected results 
and indicate that the intrinlslc dimensionality of the .A to 1A ym 
range is about four for vegetated scenes, which is what was observed 
early in the remote sensing research using twelve-band aircraft MSS 
data. 


The general conclusion on data utility is that no severe ’problems 
exist In any of the TM bands. The thermal band appears to be nu.. egls- 
tered to the right and down by four 30M pixels and this should be 
investigated further. Slight misregistrations of the 1.55 - 1.75 and 
2.08 - 2.35 bf.nds may exist. The data should be used with these prob- 
lems in mind. 
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Table 3. Calibration Accuracy Specifications for TM Br^xSs. 


Band 

Accuracy 
Specif leaf 'on 

Spec. X 
Grand Mean 

Channe 1-to -Channel 
Error Tolerance 

1 

.8) NEP 

.61 

.15 

2 

,5* ’IEP 

.16 

.04 

3 

.5% NEP 

.14 

.04 

4 

.5 % NEP 

.54 

.13 

5 

1.0* NEP 

.30 

.20 

6 

2.4* NEP 

3.42 

.85 

7 

.5° Y NETD 

MM • 

.24 •• 

.06 


** Assuming a scene mean temperature of 294° K 
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Table 3a. Arkansas TM Data. 
Stats of All Dateetors in Each Band. 


BAND 1 





Dev, of 

Dev. of 

Detector 

Mean 

Std. Dev. 

Mean 

Std. 

1 

76.35 

5.71 

-0.05 

0.02 

2 

76.29 

5.51 

-0.11 

- 0.22 

3 

76.33 

5.64 

-0.07 

-0.09 

4 

76.37 

5.77 

-0.03 

0.04 

5 

76.35 

5.78 

-0.05 

0.05 

6 

76.30 

5.65 

-0.10 

-0.08 

7 

76.35 

5.56 

-0.05 

-0.17 

8 

76.43 

5.66 

0.03 

-0.07 

9 

76.46 

5.82 

0.06 

0.09 

10 

76.48 

5.86 

0.08 

0.13 

11 

76.52 

5.93 

0.12 

0.20 

12 

76.49 

5.86 

0.09 

0.13 

13 

76.45 

5.80 

0.05 

0.07 

14 

76.46 

5.84 

0.06 

0.11 

15 

76.39 

5.56 

-0.01 

-0.08 

16 

76.37 

5.70 

-0.03 

-0.03 

Total 

1222.39 

91.74 




Grand Mean: 1222.39/16 s 76. HO 

Aver. Std. Dev.: 91.74/16 s 5.73 
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Tabla 3b. Arkansaa It* Data. 


BAND 2 


Datactor 

Maan 

Std. Day. 

Day. of 
Maan 

1 

32.95 

3.89 

0.01 

2 

32.90 

3. 78 

-0.09 

3 

32.90 

3.83 

-0.09 

9 

32.38 

3.88 

-0.06 

5 

32.90 

3.95 

-0.09 

6 

32 . 9 ; 

3.88 

•0.01 

T 

32.91 

3.83 

•0.03 

8 

32.93 

3.89 

-0.01 

9 

32.50 

3.95 

0.06 

10 

32.99 

3.99 

0.05 

11 

32.52 

3.95 

0.08 

12 

32.59 

3.97 

0.10 

13 

32.98 

3.99 

0.09 

19 

32. 9G 

3.96 

0.02 

15 

32.91 

3.86 

-0.03 

16 

32.91 

3.89 

-0.03 


Total 519. '1 62.39 


Grand Haant 519.11/16 • 32.99 
Avar. Std. Day. » 




Day. of 
Std . 

- 0.01 

-0.lt 

-0.0T 

-o.ot 

0.05 

-O.Ot 

-O.OT 

- 0.06 

0.05 

0.09 

0.05 

O.OT 

0.09 

0.06 

-0.09 

- 0.01 


62.39/16 • 3.90 
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Tab 1« 3o* Arkansas TM Data. 


BAND 3 





Dev. of 

Dev.< 

Dateotor 

Mean 

Std. Dev. 

Mean 

Std, 

1 

27.99 

6.48 

- 0.03 

- 0.10 

2 

27.95 

6.37 

- 0.07 

- 0.21 

3 

27.95 

6.46 

- 0.07 

- 0.12 

4 

27.92 

6.55 

- 0.10 

- 0.03 

5 

27.97 

6.62 

- 0.05 

0.04 

6 

28.01 

6.57 

- 0.01 

- 0.01 

7 

28.02 

6.50 

0.00 

0.08 

8 

28.03 

6.56 

0.01 

- 0.02 

9 

28.11 

6.72 

0.09 

0.14 

10 

28.13 

6.76 

0.11 

0.18 

11 

28.15 

6.77 

0.13 

0.19 

12 

28.12 

6.71 

0.10 

0.13 

13 

28.04 

6.64 

0.02 

0.06 

14 

28.07 

6.64 

0.05 

0.06 

15 

27.96 

6.50 

- 0.06 

- 0.08 

16 

27.91 

6.45 

- 0.11 

- 0.13 

Total 

448.33 

105.3 




Grand Mean: 448 . 33/16 * 28.02 

Aver . Std . Dov.i 105 . 3/16 * 6.58 
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BAND 4 


Table 3d. Arkansas TM Data. 


ector 

Mean 

Std. Dev. 

Dev. of 
Mean 

Dev.t 

Std 

1 

107.83 

21.69 

0.42 

-0.11 

2 

107.90 

21.54 

0.49 

-0.26 

3 

107.80 

21.58 

0.39 

-0.22 

4 

107.68 

21.55 

0.27 

-0.25 

5 

107.57 

21.50 

0.16 

-0.3 

6 

107.43 

21.58 

0.02 

-0.22 

7 

107.30 

21.67 

-0.11 

-0.13 

8 

107.22 

21.75 

-0.19 

-0.05 

9 

107.23 

21.93 

-0.18 

0.13 

10 

107.17 

21.99 

-0.24 

0.19 

11 

107.08 

22.01 

-0.33 

0.21 

12 

107.11 

22.06 

-0.3 

0.26 

13 

107.18 

22.15 

-0.23 

0.35 

14 

107.33 

22.08 

-0.08 

0.28 

15 

107.39 

21.98 

-0.02 

0.18 

16 

107.37 

*1.79 

-0.04 

-0.01 

Total 

1718.59 

348.55 




Grand Meant 1718.59/16 * 107.41 


Aver. Std. Dev.t 348.85/16 = 21.8 
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Table 3e. Arkansas TH Data. 


BAND 5 


Detector 

Mean 

Std. Dev. 

Dev. of 
Mean 

Dev. of 
Std. 

1 

80.02 

12.47 

0.19 

-0.21 

2 

79.95 

12.36 

0.12 

-0.32 

3 

79.94 

12.41 

0.11 

-0.27 

4 

79.91 

12.56 

0.08 

-0.12 

5 

79.91 

12.67 

0.08 

-0.01 

6 

79.87 

12.77 

0.04 

0.09 

7 

79.79 

12.83 

-0.04 

0.'.5 

8 

79.73 

12.80 

-0.1 

0.12 

9 

79.77 

12.82 

-0.06 

0.14 

10 

79.77 

12.86 

-0.06 

0.18 

11 

79.77 

12.97 

-0.06 

0.29 

12 

79.81 

12.90 

-0.02 

0.22 

13 

79.84 

12.79 

0.01 

0.11 

14 

79.86 

12.68 

0.03 

0.00 

15 

79.69 

12.57 

-0.14 

-0.11 

16 

79.68 

12.42 

-0.13 

-0.26 

Total 

1277.31 

202.88 




Grand Mean: 1277.31/16 = 79.83 


Aver. Std. Dev.: 202.88/16 = 12.68 


Table 3f. Arkansas TM Data 


BAND 6 (2.08 - 2.35u») 





Dev. of 

Dev. of 

Detector 

Mean 

Std. Dev. 

Mean 

Std. 

1 

29.37 

8.85 

0.02 

-0.24 

2 

29.32 

8.81 

-0.03 

-0.28 

3 

29.34 

8.86 

-0.01 

-0.23 

4 

29.37 

9.03 

0.02 

-0.06 

5 

29.35 

9.11 

0.00 

0.02 

6 

29.37 

9.25 

0.02 

0.16 

7 

29.34 

9.21 

-0.01 

0.12 

8 

29.39 

9.18 

0.04 

0.09 

9 

29.35 

9.28 

0.00 

0.19 

10 

29.38 

9.39 

0.03 

0.3 

11 

29.44 

9.41 

0.09 

0.32 

12 

29.45 

9.28 

0.1 

0.19 

13 

29.38 

9.17 

0.03 

0.08 

14 

29.31 

9.04 

-0.04 

-0.05 

15 

29.27 

8.85 

-0.08 

-0.24 

16 

29.23 

8.72 

-0.12 

-0.37 

Total 

469.66 

145.44 




Grand Mean: 469.66/16 s 29.35 


Aver. Std. Dev.: 


145.4*1/16 * 9.09 
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Table 4. Principal Components Results for Arkansas TM Data. 


Results Based on Bands 1, 2 , 3> 4 


EIGENVALUE 

l OF VAR. 

CUM. t 


MSE 

514.02 

89.04 

89.04 


10.96 

59.71 

10.34 

99.38 


0.62 

2.20 

0.38 

99.76 


0.24 

1.38 

0.24 

100.00 


-0.00 


EIGENVECTORS 



-0.14020 

0.59589 

-0.78599 

-0.C86.49 

-0.08969 

0.39904 

0.40832 

-0.81609 


-0.18197 

0.65220 

0.46407 

0.57110 


0.96912 

0.24560 

0.01122 

0.01920 



Results 

Based on 

Bands 1 Through 6 


EIGENVALUE 

% OF VAR. 

CUM. % 

"SE 

531.55 

64.98 

64.98 

35.02 

255.59 

31.24 

96.22 

3.78 

22.15 

2.71 

98.93 

1.07 

5.22 

0.64 

99.57 

0.43 

2.28 

0.28 

99.85 

0.15 

1.27 

0.15 

100.00 

-0.00 
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i 


EIGENVECTORS 


-0.14702 

0.22285 

0.55788 

-0.24762 0.73151 

0.14515 

-0.09673 

0.15352 

0.34346 

-0.22426 -0.54468 

0.70862 

-0.19021 

0.25512 

0.53271 

-0.05115 -0.40691 

-0.66840 

0.94779 

0.19079 

0.23642 

0.09562 -0.01416 

-0.00716 

0.00501 

0.77333 

-0.47996 

-0.40977 0.00471 

-0.06029 

-0.18579 

0.47669 

0.02749 

0.84186 0.04924 

0.16232 


Results Based 

on All Bands 


EIGENVALUE 

t OF VAR. 

CUM. t 

MSE 


525.81 

63.29 

63.29 

36.71 


253.98 

30.57 

93.87 

6.13 


24.89 

3.00 

96.86 

3.14 


18.08 

2.18 

99.04 

0.96 


4.69 

0.57 

99.60 

0.40 


2.03 

0.24 

99.85 

0.15 


1.28 

0.15 

100.00 

-0.00 


EIGENVECTORS 


-0.15027 

0.21962 

0.44266 

0.35203 

-0.27705 

-0.72814 

0.04857 

-0.09889 

0.14870 

0.28840 

0.18877 

-0.19872 

0.45909 

0.77397 

-0.19642 

0.24846 

0.44461 

0.28585 

-0.03336 

0.49591 

-0.61095 

0.93437 

0.21823 

0.10208 

0.24650 

0.08734 

0.02031 

-0.01032 

-0.00834 

0.77004 

-0.30538 

-0.38372 

-0.40432 

0.01311 

-0.05322 



-0.19541 0.47233 0.00950 0.03598 0.83933 -0.10347 0.14885 

-0.13326 0.06513 -0.64763 0.741.24 -0.08404 0.04271 -0.01591 


